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Abstract 

We introduce the Macaulay2 package BooleanGB, which computes 
a Grobner basis for Boolean polynomials using a binary representation 
rather than symbolic. We compare the runtime of several Boolean 
models from systems in biology and give an application to Sudoku. 

1 Introduction 

Buchberger's algorithm will theoretically compute a Grobner basis for any 
ideal of a multivariate polynomial ring with rational coefficients. However, 
due to memory constraints, many examples still cannot be computed in real 
time. Since all computations are done symbolically over the rational num- 
bers, two difficulties occur. The coefficients of the polynomials during the 
computation can grow very large, and also the degrees of the polynomials 
can grow very large allowing a large number of monomials in each poly- 
nomial. If an ideal consists of Boolean polynomials in the quotient ring, 
QR = ¥2[x\, . . . ,x n ]/(x\ + x%, x\ + X2, ■ ■ -x\ + x n ), we can modify Buch- 
berger's algorithm so that these two issues are avoided entirely. These ideas 
have been implemented and discussed before in [U[TU1[3]. In this paper, we 
describe a fast Grobner basis implementation for Boolean polynomials in 
C++ which is included in the Macaulay2 release version 1.4 [8]. 



2 Algorithm 



2.1 Binary Representation 

Consider a polynomial in the ring QR = F2[x\, . . . , x n ]/(xf + x\,x\ + 
X2, ■ ■ ■ x\ + x n ). Since the Boolean polynomials are multilinear, we can rep- 
resent a monomials in n variables as a binary string of length n. We store 
these strings as integers. Polynomials, then, are just lists of these integers. 
Arithmetic on the binary representation of these polynomials is very fast. 
To add two monomials, we can just concatenate the binary strings. If we 
add two identical monomials, they sum to zero. To multiply two monomials, 
we use the binary operation OR. Since xf = Xi for every i, if X{ appears in ei- 
ther monomial, then X{ will appear in the product. For the Boolean Grobner 
basis algorithm, we only need to multiply a monomial by a polynomial. To 
do this, we just multiply the monomial by each monomial in the polynomial 
using bitwise OR. Then we add the monomials by concatenation, removing 
duplicate pairs. 

The other operations that we need for Buchberger's algorithm are least 
common multiple and division of monomials. The least common multiple 
of two monomials is again just bitwise OR. In fact, the operation of least 
common multiple in this Boolean ring is just multiplication. To divide one 
monomial by another, we must first check divisibility. If a monomial m\ is 
divisible by a monomial 777,2, then any variable appearing in 777,2 must also 
appear in m\. So 7771 is divisible by 7772 if and only if bitwise m\ — 777,2 > 0. 
Once divisibility has been established, then division is binary exclusive OR. 
Note that this is nothing more than subtraction. In order to implement 
Buchberger's first criterion, we also need an operation to determine if two 
monomials are relatively prime. This can be determined by comparing the 
minimum value for each variable in the monomial. If the minimum for each 
variable is 0, then the monomials are relatively prime. Now we have enough 
operations to execute Buchberger's algorithm. 

2.2 Buchberger's Algorithm 

Given a set of Boolean polynomials, F in QR = ^[^l, • • • , x n ] / {x\ + x\, x\ + 
x<i-,---x\ + x n ), we want to compute a Grobner basis for the ideal gen- 
erated by F using Buchberger's Algorithm. The basic Buchberger Algo- 
rithm computes S— polynomials for every pair of polynomials in F. It is 
important to note that when working in a quotient ring, we need to per- 
form Buchberger's algorithm on the ideal generated by F together with the 
field polynomials. For example, consider the ideal generated by xy + z in 
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¥2[x, y, z]/(x 2 + x,y 2 + y,z 2 + z). Since there is just one polynomial in 
the ideal, Buchberger's algorithm produces the the Grobner basis for this 
ideal as xy + z. However, xz + z is also in this ideal, since xz + z = 
x(xy + z) + (xy + z). In fact, the Grobner basis for the ideal generated by 
xy + z in ¥i[x, y, z]/(x 2 + x,y 2 + y, z 2 + z) is {xy + z,yz + z,xz + z}. xz + z 
is the reduction of the S-polynomial of xy + z and x 2 + x. 

So the first difficulty in implementing Buchberger's algorithm using the 
binary representation of Boolean polynomials is how to encode these quadratic 
field polynomials. This difficulty is bypassed using the following result. 



Theorem 2.1. The S-polynomial of any multilinear Boolean polynomial and 
a quadratic field polynomial is always multilinear. 

Proof. Let / be a mutilinear Boolean polynomial with leading term lt(/). If 
lt(/) and Xi are relatively prime, then the S-polynomial of / and xf + Xi will 
reduce to by Buchberger's first criterion [5] . So we only consider quadratic 
field polynomials whose variable in the leading term is contained in ft(/). 
Supposes is in It (/). Then S(f,x 2 + Xi ) = f '+C^ 1 )(x 2 +x i ) = f-\t{f)+ Xi 
which is multilinear. □ 

Therefore, no special encoding for quadratic field polynomials is neces- 
sary. 

2.3 Implementation 

The package BooleanGB contains the algorithm gbBoolean, which is an im- 
plementation of Buchberger's algorithm for binary representations of Boolean 
polynomials over the quotient ring, QR = fl^xl, x2, . . . , x n ]/ (x 2 + x\, x\ + 
X2,x 2 l + x n ). The algorithm is implemented in C++ and is part of the 
Macaulay2 engine. gbBoolean computes a reduced Grobner basis in lexico- 
graphic order for an ideal of Boolean polynomials in QR. All computations 
are done bitwise instead of symbolically, by representing a monomial as bits 
in an integer and a polynomials as a list of monomials. On a 64 bit ma- 
chine, the algorithm works in the ring with up to 64 indeterminates. The 
algorithm uses Buchberger's first and second criteria [5] as described in [BJ 



p. 109]. Here is algorithm computing the example from Section 2.2 

11 : loadPackage "BooleanGB"; 

12 : R = ZZ/2[x,y,z] ; 
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13 : QR = R/ideal(x~2-x, y~2-y, z~2-z) ; 

14 : I = ideal (x*y+z) ; 

04 : Ideal of QR 

15 : gbBoolean I 

05 = ideal (x*y + z, y*z + z, x*z + z) 
o5 : Ideal of QR 

The code for gbBoolean is freely available with the source code distribu- 
tion of Macaulay2 [8j . 

3 Applications 

Boolean Grobner basis algorithms such as gbBoolean can be used for any 
system of Boolean polynomials. In particular, exact cover problems, satis- 
fiability problems and problems in systems biology described below are ex- 
amples. In some cases, it may be possible to turn a system of non-Boolean 
polynomials into a Boolean system. This involves dramatically increasing 
the number of variables. But in the example of Sudoku, described below, the 
time saved by the bitwise computations outweighs the increase in variables. 

3.1 Boolean Models in Systems Biology 

Logical models [1] are widely used in systems biology. They can be trans- 
lated to polynomial dynamical systems [12], in particular, logical models 
with binary variables result in Boolean polynomials. Key dynamic features 
of logical models such as fixed points correspond to the points in algebraic 
varieties generated by the polynomials describing the model. To assure that 
gbBoolean is efficient on "practical" ideals, we translate all binary logical 
models in the GINsim repository [9] to Boolean polynomial systems using 
ADAM [13] , and compute the Grobner basis of the ideal describing the fixed 
points. For all logical models in the repository, gbBoolean is faster than cur- 
rent Macaulay2 implementations, run-times are are depicted in Fig. [T] and 
table [U 
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3.2 Sudoku 



Sudoku is a popular game played on a 9 x 9 grid where the numbers 1-9 are 
filled in so that each number appears exactly once in each row column and 
3x3 block. These constraints can be represented by polynomials, with one 
variable for each cell in the grid. (See |2J and [7] for more details.) For ease 
of demonstration, we use Shidoku as an example. Shidoku is played with 
the same rules as Sudoku, but on a 4 x 4 grid with the numbers 1-4. 

For Shidoku, we use 16 variables that can each take on only the values 1- 
4. We can represent this fact, together with the constraints of the game, in a 
total of 40 polynomials. If we consider the ideal generated by these polyno- 
mials in <Q>[xi, . . . , Xis] with the lexicographical ordering, we cannot compute 
the Grobner basis for the ideal in Macaulay2 in a reasonable amount of time. 
We can, however, convert the problem to a system of Boolean polynomials 
and use gbBoolean to compute the Grobner basis. 

The Boolean system has 64 variables, representing each of the possible 
values 1-4 for each of the 16 cells. We consider the ideal of the polynomials 
in QR = F%[xl, x2, . . . , x n ]/ (x\ + x\, x\ + X2, x\ + x n ). The Boolean system 
gives us a total of 256 polynomials to represent the Shidoku constraints. 
As seen in Section [4] below, gbBoolean computes the Grobner basis for the 
Boolean ideal in just 4.1 seconds. 

4 Results 

We test gbBoolean on several random ideals and compare the run times with 
the standard Grobner basis implementation in Macaulay2. In addition, we 
also test published Boolean logical models stemming from biological systems 
and compute the basis of the ideals corresponding to key dynamic features. 
Finally, we test the algorithm on Sudoku, a problem that can be described 
by a system of Boolean polynomials. 

We compare gbBoolean to the symbolic computation in lexicographic 
order. Since the graded reverse lexicographic (gRL) monomial order usually 
yields the fastest calculation, we also compare gbBoolean to the combined 
times of computing a basis in gRL and lifting it to the quotient ring with 
lexicographic order. 

Fig. [I] shows the run-times on a 3.06 GHz Intel Core 2 Duo MacBook 
Pro. In almost all examples, gbBoolean is faster than current implemen- 
tations in Macaulay2. Cumulatively, gbBoolean is about four times faster 
( 0.27614741% ) than the Grobner basis calculation in the ring with gRL 
order, and 50 times faster ( 0.017562079% ) than the calculation in the 
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ring with lexicographic order. We do not list run-times for symbolic com- 
putations using the Sugar strategy [TT] because it is much slower than the 
"sugarless" strategy for all examples. 





QR Lex 


QR Lift 


gbBoolean 


no 


220.681 


19.0562 


4.4871 


S2 


0.002174 


0.003237 


0.00001 


S3 


0.0032 


0.004926 


0.000009 


SS2 


0.003438 


0.005541 


0.00001 


Shidoku 


470.703 


19.6797 


4.11442 


TCR 


0.007104 


0.012106 


0.000208 


THBool 


0.001946 


0.003759 


0.000009 


BoolCC 


0.002007 


0.003859 


0.000012 


erbb2 


0.002609 


0.003017 


0.00001 


yeastLi 


0.002593 


0.002013 


0.00009 



Table 1: Run Times for Random, Biological, and Shidoku Examples 




Figure 1: Benchmarks: The left graph depicts run-times for random ideals 
and an application to Sudoku, the right graph for Boolean ideals originating 
from biological systems. 



5 Conclusion 

Computing Grobner bases for Boolean polynomials bitwise rather than sym- 
bolically can speed up computation time considerably. In the example of 
Shidoku, even though converting the problem from Q to F2 increases the 
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number of variables from 16 to 64, the time to compute the Grobner basis 
decreases. The package BooleanGB in Macaulay2 has proven to be gener- 
ally faster than the standard algorithm in Macaulay2 over the quotient ring 
for Boolean polynomial ideals. Further improvements might be achieved by 
implementing the Sugar strategy or other monomial orders. 
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